Towards computability of trace distance discord 
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It is known that a reliable geometric quantifier of discord-like correlations can be built by employing the so- 
called trace distance. This is used to measure how far the state under investigation is from the closest "classical- 
quantum" one. To date, the explicit calculation of this indicator for two qubits was accomplished only for states 
such that the reduced density matrix of the measured party is maximally mixed, a class that includes Bell- 
diagonal states. Here, we first reduce the required optimization for a general two-qubit state to the minimization 
of an explicit two-variable function. Using this framework, we show next that the minimum can be analytically 
worked out in a number of relevant cases including quantum-classical and X states. This provides an explicit 
and compact expression for the trace distance discord of an arbitrary state belonging to either of these important 
classes of density matrices. 

PACS numbers: 03.65.Ud, 03.67.-a 



I. INTRODUCTION 



The issue that the quantum correlations (QCs) of a compos- 
ite state are not entirely captured by entanglement (as formerly 
believed) has recently emerged as a topical subject calling for 
the introduction of new paradigms. Despite early evidence of 
this problem was provided over a decade ago |2l an impres- 
sive burst of attention to this matter has developed only in the 
last few years [T]. In this paper, we focus on those correla- 
tions that are associated to the notion of quantum discord fSJ. 
Following the introduction of this concept, a variety of dif- 
ferent measures of QCs have been put forward (see Ref. |[T1 
for a comprehensive review). A major reason behind such a 
proliferation of QCs indicators stems from the typical diffi- 
culty in defining a reliable measure that is easily computable. 
No general closed formula of quantum discord, for instance, 
is known (with strong indications that this is an unsolvable 
problem lU) even for a pair of two-dimensional systems or 
"qubits" pll, namely the simplest composite quantum system. 
Unfortunately, the demand for computability typically comes 
at the cost of ending up with quantities that fail to be bona fide 
measures. In this respect, the most paradigmatic instance is 
embodied by the so called geometric discord (GD) Q, which 
while being effortlessly computable (and in some cases able 
to provide useful information) suffers from a major drawback. 
It can indeed grow under local operations on the unmeasured 
party |'6l, an effect which a bona fide indicator (e.g. quantum 
discord) is required not to exhibit. Following an approach fre- 
quently adopted for other QCs measures, the one-sided GD 
is defined as the distance between the state under study and 
the set of classical-quantum states. The latter class features 
zero quantum discord with respect to the measured party, say 
subsystem A, which entails the existence of at least one set 
of local projective measurements on A leaving the state un- 
perturbed |T "71). While the above definition in terms of a 
distance is clear and intuitive, it requires the use of a metric 
in the Hilbert space. The GD employes the Hilbert-Schmidt 
distance, which is defined in terms of the Schatten 2-norm. 



Such a distance is well-known not to fulfill the property of be- 
ing contractive under trace-preserving quantum channels |8|, 
which is indeed the reason behind the aforementioned draw- 
back of GD L9J. This naturally leads to a redefinition of the 
GD in terms of a metric that obeys the contractivity property. 
One such metric is the trace distance f3i|, which employs the 
Schatten 1-norm (or trace norm for brevity). In the remainder 
of this paper, we refer to the QCs geometric measure resulting 
from this specific choice as trace distance discord (TDD). 

While investigations are still in the early stages lfT0lfT4l . 
TDD appears to enjoy attractive features. Besides the dis- 
cussed contractivity property, the trace distance is invariant 
under unitary transformations. More importantly, it is in one- 
to-one correspondence with one-shot state distinguishability 
fT5l|, i.e., the maximum probability to distinguish between 
two states through a single measurement. This operational 
interpretation provides evidence that the trace distance works 
as an accurate "meter" in the space of quantum states. An- 
other appealing advantage of TDD lies in its connection with 
entanglement. Recently, indeed, it was suggested to define 
the full amount of discord-like correlations in a system S as 
the minimum entanglement between S and the measurement 
apparatus created in a local measurement (see ifTSU TBl and 
references therein). This way, a given entanglement measure 
|fT9l identifies a corresponding QCs indicator. Remarkably, it 
turns out that the latter always exceeds the entanglement be- 
tween the subparts of S when this is quantified via the same 
entanglement measure. This rigorously formalizes the idea 
that a composite state can feature QCs that cannot be ascribed 
to entanglement. In this framework, it can be shown lfT2ll that 
the entanglement counterpart of TDD is negativity ['20 1, the 
latter being a well-known - in general easily computable - 
entanglement monotone 11211 . 

In spite of all such interesting features, the easiness of com- 
putation of TDD in actual problems is yet to be assessed. To 
date, the only class of states for which a closed analytical ex- 
pression has been worked out are the Bell-diagonal (BD) two- 
qubit states, or more generally states that appear maximally 



mixed to the measured party irT2l[T3l . While the proof of this 
formula is non-trivial ||12j. this does not clarify whether or 
not, besides its reliability, TDD brings about computability 
advantages as well. Owing to the high symmetry and reduced 
number of parameters of BD states, indeed, most if not all of 
the bona fide QCs measures proposed so far can be analyti- 
cally calculated for this specific class f23]. 

In this paper, we take a step forward and set up the prob- 
lem of the actual computation of two-qubit TDD on a new 
basis. We first develop a theoretical framework that reduces 
this task to the equivalent minimization of a two-variable ex- 
plicit function, which parametrically depends on the Bloch 
vectors of the marginals and the singular values of the cor- 
relation matrix. Next, after re-deriving the value of TDD for 
a class of density matrices that includes BD states, we discuss 
two further relevant cases in which the minimization problem 
can be analytically solved. One is the case where the cor- 
relation matrix has one non-zero singular eigenvalue, a sub- 
set of which is given by the quantum-classical states (unlike 
classical-quantum states these feature non-classical correla- 
tions with respect to party A). The other case is given by 
the family of X states |22|, which include BD states as spe- 
cial cases. While these are arguably among the most stud- 
ied classes of two-qubit density matrices |1 1, the calculation 
of their QCs through bona fide measures is in general a de- 
manding task. To the best of our knowledge, in particular, no 
closed expression for an arbitrary quantum-classical state is 
known to date with the exception of Ref. Il24l where however 
an ad hoc measure exclusively devised for this specific class 
of states was presented. In the general case, indeed, one such 
state depends on four independent parameters and, moreover, 
features quite low symmetry. In Refs. f24l and ESJI . for in- 
stance, closed expressions for a fidelity-based measure | |26) 
and the quantum discord, respectively, could be worked out 
only for high-symmetry two-parameter subsets of this family. 

Even more involved is the calculation of QCs in the case 
of X states, a class which depends on five independent pa- 
rameters. Regarding quantum discord, an algorithm has been 
put forward by Ali et al. Il27ll . Later, however, some coun- 
terexamples of X states for which such algorithm fails were 
highlighted |28| (see also Ref. 1 1 1). 

The present paper is organized as follows. In Sectionlll] we 
present our method for tackling and simplifying the calcula- 
tion of TDD for an arbitrary two-qubit state. This is demon- 
strably reduced to the minimization of an explicit two-variable 
function. In Section III we apply the theory to the case of 



Bell states and that of density matrices having correlation ma- 
trix with uniform spectrum. In Section |IV] we show that the 
minimum can be analytically found in a closed form whenever 
the correlation matrix of the composite state features only one 
non-zero singular value. As an application of this finding, in 
Subsection |IV A| we compute the TDD of the most general 
quantum-classical state. As a further case where the mini- 
mization in Sectionlnlcan be performed explicitly, in Section 
M we tackle the important class of X states and work out the 
TDD for an arbitrary element of this. We finally draw our con- 



clusions in Section VI A few technical details are presented 
in the Appendix. 



II. ONE-SIDED TDD FOR TWO-QUBIT STATES: 
GENERAL CASE 

The one-sided TDD D^^^ipAs) from A to B of a bipartite 
quantum state pab is defined as the minimal (trace norm) dis- 
tance between such state and the set CQ of classical-quantum 
density matrices which exhibit zero quantum discord with re- 
spect to local measurements on A, i.e. states which admit an 
unravelling of the form 



Pab^ = X \"J^A{aj\'^eB(j) 



(1) 



with 10";)^ being orthonormal vectors of A and QBij) being 
positive (non necessarily normalized) operators of B. Specif- 
ically, if ||0||i - Tr[ V©'©] denotes the trace norm (or Schat- 
ten 1-norm) of a generic operator © then 



1 



"D^^Hpab) = T min \\pab 



W«\ 



-p'abWi 



(2) 



the 1/2 factor ensuring that D'^^Hpab) takes values between 
and 1 [analogous definition applies for the one-sided TDD 
from B to A, D^'^Kpab)]- The quantity in Eq. ^ fulfills 
several requirements which make it fit for describing non- 
classical correlations of the discord type [12]. In particular, 
from the properties of the trace distance Ol it follows that 
D^^^(Pab) 

i) is zero if and only if pab is one of the classical-quantum 
density matrices ([TJ; 

ii) is invariant under the action of an arbitrary unitary op- 
eration Ua ® Vb that acts locally on A and B, i.e. 



= Tl(^) 



n^-^'iPAB) = n'^'iUA ® Vb pab f/; ® Vfi) ; 



(3) 



Hi) is monotonically decreasing under completely positive 
and trace preserving (CPT) maps on B; 

iv) is an entanglement monotone when pab is pure. 

Furthermore, in the special case in which A is a qubit Eq. (J2]) 
can be expressed as lfT2ll 

S'^Vab) = r min \\pAB - (Ua ® ]1b)(pab)IIi , (4) 

Z [Ha] 

where now the minimization is performed with respect to all 
possible completely depolarizing channel 11^ on A associated 
with projective measurements over an orthonormal basis, i.e. 



Ua(-)^Pa-Pa + Qa-Q 



(5) 



withPyi = |*I')yi(*I'| and 2a - 1a-^a being rank-one projectors 
(I*!*) is a generic one-qubit /?Mre state). 

In what follows, we will focus on the case where both A 
and B are qubits. Accordingly, we parametrize the state pab 



in terms of the Pauli matrices (?a(B) = {o-a(B)uo-a(B)2, <7'a(B)3) = 

{o'A(B)x, 0-A(B)y, 0"a(B)z), 1-6. 



PAB = -(iA<S>iB + XA-0-A<8>tB+tA<8>XB-O-B+2^ ^ ija-Ai®CrBj) 

(6) 



U=l 



where 



XA{B) - Tr[pABCr^(B)] 



(7) 



is the Bloch vector corresponding to the reduced density ma- 
trix Pa(B) describing the state of A{B), while F is the 3x3 real 
correlation matrix given by 



Tij = Tx\pAB{o-Ai ® a-Bj)] 



(8) 



Similarly, without loss of generality, we express the orthogo- 
nal projectors Pa and Qa of Eq. (|5]l as 



1 



-(1a +e- &a) 



Qa 



1 



-{lA-e-&A) (9) 



with e being the 3-dimensional (real) unit vector associated 
with the pure state |*I')a in the Bloch sphere. Using this and 
observing that nA(lA) - 1a, and liAiv ■ o^a) -(e ■ v)(e ■ o-a), 
Eq. (pl can be arranged as 



0^^\pab) 



1 



■min||M(e) 



(10) 



where the minimization is performed over the unit vector e 
and M{e) is a 4 x 4 matrix which admits the representation 

M(e) = [{xa - (e ■ XA)e) ■ (?a] ® Ib 

+ 2_j ^U ^^' ~ ^'^) -d-A^o-Bj ■ (11) 

ij 

Here, i, is the /th Cartesian unit vector and e, = xi-e the /th 
component of e (note that ctaj - i, ■ o-a). The second term 
in Eq. ( [TT| can be further simplified by transforming F into a 
diagonal form via its singular value decomposition I.29J . More 
precisely, exploiting the fact that F is real we can express it as 



F = 0^ 



n 





1 





72 











73 _ 



Q, 



(12) 



where O and Q are real orthogonal matrices of SO(3) while 
{ji] are real (not necessarily non-negative) quantities whose 
moduli correspond to the singular eigenvalues of F Il30l . We 
can then define the two set of vectors 



7=1 ;=i 



(13) 



for k = 1,2,3. As (9, Q e SO(3), by construction {w^} is an 
orthonormal (right-hand oriented) set of real vectors and so 
is {vic} (each is indeed a rotation of the Cartesian unit vectors 



{xj]). Using the above, we can arrange Eq. ( 1 1 



I as 



M(e) = (xAi. ■ o-a) » 1b + y 7k (wt± ■ (?a) ® {vk ■ <?b) , (14) 



jt=i 



where for compactness of notation we introduced the vectors 

XAj_^XA-(e ■ XA)e , Wk±^Wk-(e-Wk)e (15) 

to represent the orthogonal component of xa and Wk with re- 
spect to e. 

Note that {iik ■ (?b) describes the transformed set of Pauli 
matrices under a local rotation on B. This set clearly fulfills 
all the properties of Pauli matrices as well. One can there- 
fore redefine the B's Pauli matrices as {vk ■ <?b) — > crBk, which 
amounts to applying a local unitary on B. Let then M'(e) be 
the transformed operator obtained from M(e) under such rota- 
tion, i.e. 



M' 



'(e) = (xa± ■ (?a) ® 1b + 2_] 7* ("^k^ ■ ^a) ® cTBk ■ (16) 



k=\ 



Since the trace norm is invariant under any local unitary we 
have 



||M(g)||i=||M'(g) 



(17) 



in line with the invariance property ii) of D^'^Hpab) [indeed 



M'{e) is the operator (111 associated to the state p^^ obtained 



from pab via a local unitary rotation associated to the trans- 
formation [vk ■ d'B)^' cTBk}- The trace norm of M'(e) can now 
computed by diagonalizing the operator M'{ej^ M'ie). For 
this purpose, we recall that given two arbitrary vectors [x, y], 
the Pauli matrices fulfill the following commutation and anti- 
commutation relations 

[x-d-A,y-^A\ ^ '^i(x^y)■^■A, (18) 

[x- o-A , y-d-A] ^ 2{x-y) (19) 

as well as the identities 0-^10-^2 = icTAi, ctaictai - -icta^ and 
the analogous identities obtained through cyclic permutations 
(in the above expression "A" indicates the cross product). Us- 
ing these, we straightforwardly end up with 

M'(e)"'"M'(e) = (g + x^ J I^b + A + 2 l^ ®;? ■ <?b, (20) 

where jca± = \xaA (throughout, x - \x\ for any vector x), ;t^ is a 
tridimensional real vector of components 



Xk = Jk Wk± ■ Xa± , 

while 2 is a positive quantity defined as 

3 



e = 2 ^*' 1*^*^1' 



(21) 



(22) 



k=\ 



and finally A is the operator 

A = 2_] 'yj'y^ [('^'t-l- ■ ^A)(Wk± ■ (?a) ® O-BjO-Bk] 



j*k 



- -27172 (Wii A W2A.) ■ (?A ® O-B-i 
-'^.yiji (W2j. a W3j_) ■ &A ® a-B\ 
-27371 (W3j_ a Wij_) ■ &A ® 0-B2 ■ 



(23) 



This expression can be simplified by observing that since the 
Wk±'s are vectors orthogonal to e [see Eq. ( 15 1] their mutual 
cross products must be collinear with the latter Indeed, intro- 
ducing the spherical coordinates {0, 0) which specify e in the 
reference frame defined by {w^}, we have 



(wix A vv2j_) = (m>3 ■ e) e — cos 6 e , 
(w2j_ a vv3x) = (wi ■ e) e — sin 6 cos (p e , 
(w3j^A vvij_) = (m>2 ■ e) e — sin sin e . 



(24) 



Substituting these identities in Eq. ( [23] ), the operator A can 
remarkably be arranged in terms of a simple tensor product as 

A = -2(e-(?^)®(|-<?B), (25) 

where g is the vector 

g- (7273 sin cos 0,7371 sin sin 0,71 72 cos 0) , (26) 

which is orthogonal to;^^ 131]. Next, observe that the operator 
e ■ (?A of Eq. (25 1 is Hermitian with eigenvalues 1 and -1. 
Therefore, if {10)^,11)^1 are its eigenvectors we can write e ■ 
(?a = |0)a<0|-|1)a<1|. Plugging this and 1^ = |0)a<0| + |1)a<1| 



into Eq. ( 20 1 this can be arranged as 



+ 2 [\0U0\ ® 0?-^-(?B +|1)^<1| ® (x+8)-(?b] , 

which can now be put in diagonal form. Indeed, due to the 
aforementioned spectrum of x ■ <?, it has eigenvalues A - 
Q+x^^+2 yjx^+g^^ sach twofold degenerate Ii31il . Therefore, 



through Eq. ( 17 1 we end up with 



\\M{e 



li =2 1 A/fl + V^ + A/fl - V^ ) , (27) 



where 



a = a{e) — Q+Xj^^ — Q+Xj^ — (xA-e) , 



b = b{e)^A[x'+g')- 



(28) 
(29) 



Note that Q, xaj_, x ^"d g are all functions of e [cf. Eqs. (21 1, 



(22i and(26i]. As||M(e)||i is a positive-definite function, find- 
ing its minimum is equivalent to searching for the minimum 
of its square ||M(e)||p Thereby 



mm||M(g)||i= ^mm||M(e)||2 = 2y2 

where the function h{e) is defined as 

h{e) - a{e) + ^Ja^ie) - b(e) . 



min h(e) 



(31) 



In conclusion, in the light of Eqs. (10 1, ( 27 1 and ( 30 1 



1 



^'^"Vas) = jmin 

1 



^Ja+ yfb+^a-yfb 



min/z 



(32) 



wz [73] 




W2 [72] 



FIG. 1 : (Color online) Schematics of the minimization procedure for 
calculating the trace distance discord D'^*(p^b) of a two-qubit state 
pAB- The reference frame in which this is carried out is defined by 
the orthonormal set of three vectors {wu}, where each Wk is associ- 
ated with a real singular eigenvalue 7^ of the correlation matrix [see 
Eqs. ([8}, l |12[ l and l|T3}]. This frame identifies a representation for the 
local Bloch vector x^ [defined in Eq. (171]. All these quantities are 
drawn using solid black lines to highlight that, for a given density 
matrix Pab, they are fixed. Instead, the unit vector e (red line) rep- 
resents the direction along which a projective measurement on A is 
performed. In the optimization procedure, e is varied until function 
h in Eq. (|31^ reaches its global minimum according to Eq. \'i2). 



We have thus expressed our trace-norm-based measure of QCs 
of an arbitrary state p^B as the minimum of an explicit function 
of the two angles [0, cp} (O<0<n, 0< 



I < 2;7r). Equation ( 32 1 
is the first main finding of this paper. For clarity, all quantities 
involved in the minimization problem under investigation are 
pictorially represented in Fig.fT] 



III. BELL DIAGONAL STATES AND STATES WITH 
HOMOGENEOUS SINGULAR VALUES 



The optimization in Eq. ( 32 1 simplifies when the state pos- 



sesses certain symmetries. In particular, by ordering the sin- 
gular eigenvalues of T as (this convention is adopted only in 
the present section) 



I71I > M > IrsI , 



(30) one can show that 



D'^'iPAB) 



\n\ 



(33) 



(34) 



at least for two classes of states p^s, which we label as 'A' and 
'B', respectively. These are defined as 

class A : arbitrary {7^} but ^a = 

class B : arbitrary xa but {7^) with the equal moduli, i.e. 

|7,|=7 V^= 1,2,3. (35) 

We develop the proof in the following two subsections. 



A. Bell diagonal states 

States of class A, which include Bell diagonal states, are 
characterized by the property that the reduced density matrix 



of subsystem A is maximally mixed. For these, Eq. ( 32 1 was 
proven in Refs. lfT2l[T3l using an independent approach. Here, 
we present an alternative (possibly simpler) derivation based 



on Eq. ( 32 1. We point out that these states form a special sub- 



set of X states, which we will study in full detail in SectionlV] 
Here, our goal is indeed to present a straightforward applica- 
tion of our method for calculating the TDD developed in the 
previous Section. 

To begin with, we observe that if xa - then the vector 
X in Eq. (21 1 vanishes, i.e. ;^ = 0, while the function a in 



Eq. (28 1 coincides with Q in Eq. (22i. Thereby, the function 



h in Eq. ( 3 1 1, which we have to minimize over e according to 
Eq. (p2l), becomes 



h^Q + ypH with H^Q^ - Ag^ 



(36) 



Expressing now Q in terms of 6 and 4> and due to the ordering 



in Eq. ( 33 1, it turns out that 



Q{e, 4>) = yf ( 1 - cos V sin^ 0) -H 72( 1 - sin^ sin^ 9) (37) 



+ yl(\ -cos^6')> Q{e = Til2, 



0)^yl+yl 



namely Q reaches its minimum value for 6 - 7r/2 and (p - 0, 
i.e. when e points toward vvi. The same property holds for the 
function H. Indeed one has 



H(e, 4>) = A{0) sin'* (f> + 8(6) sin^ (f> + C(0) 
> H(6^7T/2,4>^0)^(yl-yj)\ 



(38) 



where we used 



A(e) = sin^ (y\ - ylf > A(0) = , 
B(9) = 2 sm^0(y]- y\)\{y\- y\) 

H- cos^ {y\ - y\)\ > 8(0) = , 
C(0) = (yl-ylf+2cos^0(y'2+yl)(yl 

.,2\2 



rl) 



2x2 



+ cos" 0(yi - yir > C(0 = n/l) > (y'^ - y^) 



Replacing Eq. (|37|) and ^Ej into Eq. (|36jl enatils h(0, (p) > 
h{0 = 71 12,4) = 0) > 2\y2\, which through Eq. (|32]) yields 
Eq. (ig. 



B. States with homogeneous lyt I 's 

Class B (see definition given above) includes, for instance, 
mixtures of the formp^^ - pg^ ® Ib/2 + (1 - p)\^-)ab{^-\ 
where p e [0, 1], Qa is an arbitrary state of A, and |*I'ab) is 
the singlet state ^Vab) = (|01)ab-|10)ab)/V2 [from now on, 
{|0)a(B), |1)a(B)) denotes an orthonormal basis for A (B)]. In 
this case, yt = (1 - p) for all k while xa - psA with sa the 



Bloch vector of ^i^: therefore according to Eq. (34i this state 
has a value for TDD given by (1 - p)/2. 



To derive Eq. ( 34 1, we introduce the diagonal matrix T = 
diag(fii, f22, ^33) formed by the coefficients fn, ^22, ^33 defined 
by the identities 



yj = 'jj y ■ 



(39) 



[it is clear from (35 1 that tjj can only take values ±1]. Under 



this condition, from Eqs. (21 1, (22 1 and (26 1 it then follows 

.,2 



Q = 2y' 



^y^T 



l^f =/ 



X ^ yTxA,^^ iri' = r'l-4,j' , (40) 

where ^ takes value either 1 or - 1 depending on the explicit 



form of the mapping ( 39 1. Replacing this into Eqs. ( 28 1, ( 29 1 



and ( 3 1 1 we end up with 



h^2y^ + 2\^A,J 



(41) 



which depends upon e through |iti,±|^ only. The minimum is 
then achieved when |xa,±I vanishes, which clearly occurs by 
taking e along the direction of xa [recall Eq. (15 1]. Thus 



min/i = 2y 



(42) 



which when replaced into Eq. ( 32 1 gives Eq. ( [34) i, as antici- 
pated. 



IV. CORRELATION MATRIX WITH A SINGLE 
NON-ZERO SINGULAR EIGENVALUE 

This class of states is important since quantum-classical 
states fall within it, as we show later. It is defined by [see 



Eq. (12 1] 72 = 73 = while yi=y and xa are arbitrary (the only 
constraint is that the resulting pab must be a properly defined 
density matrix). We show below that the TDD of one such 
state is given by 



D'-^'ipAB) 



Ifl A Xa\ 



mm 



1 



Ifl ±^a\ 



(43) 



where yi = lyilvvi, vvi being the first element of the set {wk} 



defined in Eq. (13 1. Eq. (43 1 is another main finding of this 
work. 

To begin with, we observe that due to 72 = 73 =0 we are 
free to choose the direction of the Cartesian axes vv2 and w^ 
(vi>2 ±1^3) on the plane orthogonal to vvi. We thus take vv2 as 
lying on the plane formed by vi>i and xa- Hence we can write 
Xa - xaiWi + XA2W2, where xai and xa2 are the components 
of iA in reference frame defined by {wk}- Accordingly, 



xai—xa-wi—xa cos a , xa2 — xa-W2 — xa sm a 



(44) 



with a being the angle between xa and vvi while xa 



4' 



A\ 



■X^ 



With the help of Eqs. (21 1, ( 



present case a and b [cf. Eqs. (28 1 aM 



m\ 



22 1 and (|26|), in the 
read 



a^y^ +x\-{y'-e\He ■ xa)\ , b^4y^[xAi-eiie ■ xa)]\ (45) 



where ei = e ■ m>i. Observe then that we can write 



a + ^/b -(y + icAif + i^ 



A2 



[ye I ± (e ■ Xa)] 



(46) 



It turns out that both a+ yfb and a- 'sfb decrease when the com- 
ponent of e on the plane formed by m>i and xaA-S-, the wi -W2 
plane, grows. To see this, we decompose e as e = s + Bj_, where 
e = eiM>i + e2W2 is the component of e on the W1-W2 plane, 
while s^ = e^Wi the one orthogonal to it. With these defini- 
tions, in Eq. (46 1 we can evidently replace e with e (we remind 
that XAi - 0). Now, it should be evident that the last term of 



Eq. ( 46 1 can be written as - [yei + (e ■ x^)] =-|ep[/±(<^, a)]^, 
where f±{<p, a) = y cos (p+XA cos{^-a) is a function of (p (i.e., 
the azimuthal angle of e) and the aforementioned a. Clearly, 
for given <f> the minimum of a + V^ is achieved when |e| is 
maximum, i.e., for s= e or equivalently 6 = n/l. Thus, due to 



Eq. ( 27 1, in Eq. ( 32 1 we can safely restrict the minimization 
over e - (6, (p) to the set e - (njl, (p). To summarize, we need 
to calculate 



mm 



\\M(e) 



-2^, ^{y + riXAif+xl^-f,{cl>,a)\ (47) 



';=± 



Through few straightforward steps (see Appendix A), ||M(e)||i 
can be arranged as (we henceforth omit to specificy 9-nl2) 



||M(e)||i=2^ |xa s\n{(p - a) -1-777 sin i 



(48) 



rj=± 



Exploiting the positiveness of ||M(e)||i and the identity (ly-i-zl-i- 
\y-z[f'=4-max{y^,z^}, where y and z are any two real numbers, 
Eq. (|48|) can be converted into 



||M(e)||i=4max{|xAsin(0-Q')|,|ysin<;*|) (49) 

-4-Jx^ + y^ max {\ sin /3 sm((f> - a)\,\ cos /3 sin cf)\] , 
where the angle /3 is defined through the identity 

i/3^\xA\/ylxl+y^. 



sm/ 



(50) 



Replacing ||M(e)|| so obtained into Eq. (lOi we can then ex- 
press the one-sided TDD of our state pab in terms of the fol- 
lowing min-max problem, 

D^-^(PAS)-^ 



X min max{| sin/?sin(0-Q')|, |coSj6sin((;A)|). (51) 

(^G[0,2;r] 

An analytic solution is obtained by observing that the (f>- 
dependent functions /i(^) - |sinySsin(0 - a)\ and /2(0) - 
|cosy6sin(0)| have the same period n and that in the domain 
0€ [0, tt] exhibit the two crossing points <pc+ and 0c_ given by 



cot(0c±) = cot a : 



cot/3 



smff 



(52) 



By construction, the function Eq. (49i reaches is minimum 
either in (pc+ or in (pc-. Therefore, 



n'^'iPAB) 



V^ 



+y- 



2 
\7XA2\ 



min{|coSjSsin((;*c+ 

1 



nun 



^(■^±1-4 1)^1-^2 



|cosy6sin(0(;-)l) 
(53) 



where the latter identity have been obtained through simple 
algebraic manipulations. To arrange this formula in a form 



independent of the reference frame, we make use of Eqs. ( 44 1 



and (50 1. This finally yields Eq. (43 1 



A. Quantum-classical states 

The result of the previous section can be exploited to pro- 
vide an analytical closed formula of D^^^(pab) for the well- 
known class of quantum-classical states. One such state reads 

Pab ^PPoa® |0>b<0| + {I - p)piA ® |1>b<1|, (54) 

where po(i) is an arbitrary single-qubit state with associated 
Bloch vector iod), i.e, pod) = {t+S{x\)-(f)l2. The state in 



Eq. (54 1 represents a paradigmatic example of a separable 
state which is still able to feature A — » B quantum correla- 
tions. On the other hand, note that the quantum discord in the 
opposite direction, B— »A, is zero by construction. 

One can assume without loss of generality that sq - (0, 0, ^o) 
and s\ - {s\ sin ip, 0, s\ cos ip) with < ^ < tt, i.e., the Z-axis of 
the Bloch sphere is taken along the direction of sq while the 
F-axis lies orthogonal to the plane containing both s\^ and s\ . 
Vector xa and matrix F are calculated as 

Xa - iX^-p)s\sinip,Q,psQ + {\-p)siCosip) , (55) 



r = 



' (/:i-l)iisin^ 


pso + {\- p)s \Cos ^p , 



(56) 



r has only one singular eigenvalue since its singular value de- 
composition yields 72 =73 =0 and 



lril=r= ^|p'^sl+(p-\)sl{{p-\)sx+2psQC0s^p\. (57) 
Such states therefore fall exactly in the case studied in the pre- 



vious section. To apply Eq. (43 1, though, we need to calculate 



the unit vectors Wk- From the matrix Eq. (56 1, they are calcu- 
lated as 

w\ - —({p-\)sisinip,Q,pS() + {p-\)siCOSip\ (58) 

At ^ ^ 



1 l{l-p)siCOtip-pSQCSCip 

wi = T- ; 77 ,0,1 

A2\ (p-l)si 

ws = (0,1,0), 



(59) 
(60) 



where Ai 2 are normalization coefficients. In particular, it 
turns that Ai coincides with y in Eq. ( [57] i, i.e. Ai = y. Hence, 
the vector yi -ywi in Eq. (43 1 is given by 



n 



- ((p-l)si sin<p,0,pso+(p-l)si cos(p) . (61) 



This, together with Eq. (55 1, yields the identities 



IfAAfil = 2p(l- p) sosi\sin(p\ , 

\yi + xa\ = 2p sq , 

Ifi -xa\ = 2(1 - p) si . 



Replacing these into Eq. (43 i, we end up with 



&^'(Pab) 



I sin (/3 1 



min{pio,(l - P)s\ 



(62) 



which represents the TDD of the most general quantum- 
classical state [Eq. (54i]. This has a very clear interpre- 
tation once one notices that the maximum value of D*^^ 
is 1/4 and is taken for: so - s\ - I, p - \I2 and ip - 
7t/2. For these parameters, Eq. (|54]i reduces to pab - 



1/2(|0>^<0|®|0>^0| + |+U+|®|1>^1|) (where |±> = (|0> + 
|1))/V2), which is a popular instance of separable quantum- 
correlated state. The qualitative behavior of £) for sq = ^i and 
;? = 1 /2 is fully in line with that of the quantum discord |'25| 
and for io = ■si = 1 with that of the fidehty-based measure 
analyzed in Ref. Il24l . 



subject to the constraints Yjt=i Pa - 1' P11P44 ^ Ipi4p 

> IP23P 



and 



P22P33 ^ IP23r- Here, we have referred to the computational 
basis {|00>AB, \01)ab, |10)as, |11)ab}- Without loss of general- 
ity, off-diagonal entries P32 and p4i can be taken as positive, 
i.e., p32>0 andp4i >0 |32|. It is straightforward to check that 
for such states xa,i =xa_2=0, Xy(3 =2(piiH-p22)-l (that is, ^^ 
lies along the i3-axis) while the correlation matrix has already 
the diagonal form since F^diag {71,72, 73) with 



7l=2(p32-l-p4l), 72 = 2(p32-p4l), 73 = 1 -2(p22+P33)- (64) 



V. X STATES 



A two-qubit X state has the X-shaped matrix form 



Hence, in the present case Wk - i/t for A: = 1 , 2, 3 [see Eq. ( 1 3 1] 



We therefore have to deal with the four parameters XAi and 
(7^;). Note that the only hierarchical relation which always 



holds is I7i|>l72l [seeEq. (64i]. 



Pab 



Pii p^i 

P22 PI2 

P32 P33 

P41 P44 j 



In what follows, we will prove that the TDD of state 



(63) Eq. ( 63 1 is given by 



if7?-73+43<0 



if 7^-72+43 >0 



©*^Vab)=^, 



fif|73l>l7ll r)'^^HpAB)-- 



lyil 



ifi73i<i7ii D<-v..)=0(7^r3^+43)W'^lw^^®[-(^^^^^-)l^' 



(65) 



where we have used the Heaviside step function &(x) [we 
adopt the standard convention 0(0) = 1/2]. It can be checked 
(see Appendix B) that for Bell-diagonal states the above ex- 
pression reproduces the result of Section III i.e., the TDD 



is half the intermediate value among Wykl) [we stress that 
here the labelling of the 7/:'s does not imply the ordering in 
Eq. Q]. 



Equation ( 65 1 can also be written in the compact form 



&^>(Pab) = 



1 I y] max{yl,yl + xl^]-yl min{72,72) 
2"^ max{yl,yl + xl^}-min{yl,yl}+yl-yl 



(66) 



showing that for the X-states the discord is only a function of 
the following three parameters: I71I, I73I and 



^yl + x2 



A3- 



To begin with, Eqs. (28 1 and (29i imply that the (6. 



dependent functions a and b entering Eq. (31 1 [recall that 



(6, (p) specify e] depend only on /v s sin 6 and v = sin (p as 



a-aQ+axfj., b = b[)+bi fi+b2fJ. , 



(67) 



where (a,) and {/?,) are the following linear functions of v 

ao= 7?+72, ai^iyl+xl^-y]) + iy^i-yl)v, (68) 

Z7o=47?7i b2^4xl,[(yl-yi) + (y^i-yl)v], (69) 
^1=4 [jlyl+y2^(xl,-yl) + (y',-yl)iyl-xl,)y] . (70) 

Clearly, a(fi, v) and b(fi, v) are defined in the square S defined 
by >S = {//,v: 0<//<l,0<v<l) [and so is /z = a + ^a^-b. 



see Eq. (31 1]. The partial derivative of h with respect to v, 
dvh, can be arranged as dvh - (2hdva-dyb)/ (2 y/a^-b) (an 
analogous formula holds for cJ^/i). Now, due to Eqs. ([67|)-([70|) 
dya-(y^-yl)li and, notably, dyb -4[y^+{^-l)x\^]dya~when 
these are replaced in dyh we thus end up with 



h-2[yl+xl^{t,-l)] 

Oyh - -^=^ OyO. 

Vfl2 -b 



(71) 



As witnessed by the denominator of this equation, we observe 
that function h is in general non-differentiable at points such 
that a^ - b, owing to the square root ^d^-b appearing in its 
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definition, Eq. (31 1. One then has to investigate these points 
carefiilly, as they may potentially yield extremal values of h 
that would not be found by simply imposing d^h - dyh - 0. 

As a key step in our reasoning, we first demonstrate that a 
minimum ofh cannot occur in the interior ofS. Afterward, we 
minimize function h on the boundary of S, which will even- 



tually lead to formula (65 i 



A. Proof that minimum points cannot lie in tlie interior of S 

We first address minimum points at which h is differen- 
tiable, i.e., that fulfill a^ +b entailing the existence of partial 
derivatives for h. A necessary condition for h to take a mini- 



mum on these points is then dyh - 0. Based on Eq. (71 1, this 
can happen when either /! = /!o=2[73H-x^2(jU - 1)] ox dyQ-Q. 

In the latter case, as discussed above, dyU - iy^-Jo)fJ., which 
vanishes for // = (that is on the boundary of S) or |yi | - 
|y2l- Using Eq. (31 1 and Eqs. (67i through (70i it is easy to 



calculate that when |yi| = lyal, depending on the sign of yj" 
yl+x^y either h = liy^+x^^ij.) or h = 2[y\+{yl^-y^^] . Thereby, 
in the case \y\ | = 172! the minima of h must fall on the boundary 
of>S. 

Let us now analyze the situation where h-h^- 2[y^+x^^(ju- 



1)], which would also yield dyh = [cf. Eq. (71 1]. As h-a + 
'^a^-b, a necessary condition for this to occur is clearly (/zq- 
a)^ — a^-b. With the help of Eqs. ( 67 1-([70|, this identity can 
be explicitly written as 4(l-yu) \y\-y\+x\^\y\-y\+x^^.A — Q>. 
This is fulfilled if at least one of the following identities holds: 
(i) // = 1 , (ii) y\-yl + x\^ = 0, (iii) yl-yl+xl^^O. Case (i) 
clearly corresponds to a point on the boundary of S. In case 



(ii), using that yi^ ■ 



■y\+x^ 



A3 



<0 (due toy^ <y\ =73-^^3) 



we end up with h - 2{(y\ — x^^ + x^^-^^}. In case (iii), using 

) we have that 



that yi -73 + Jc^3 > (due to y\>yi 



h^2{y} + [y^ 



-y\- 



•^A3 



3 7i +(7f ~T3+-^^3)^]A')' whose minimum occurs 



for yu = V = or ju = 1 and v = (depending on the sign of 
y\-y'\)- Hence, even in cases (ii) and (iii), the minima of h 
fall on the boundary of S. The above shows that no minima 
points at which h is differentiable can lie in the interior S. 

Let us now address singular points, i.e., those at which h 
is non-diff'erentiable and hence minimization criteria based on 
partial derivatives do not apply. These points (see above dis- 
cussion) are the zeros of the function f - a^-b. Our aim is 
proving that even such points, if existing, lie on the bound- 
ary of S. Firstly, note that / > (we recall that a^ >b al- 
ways holds, see Section 11). This means that a zero of / is 
also a minimum point for /. From Eqs. ([67|-([70li it is evi- 
dent that /(/i, v) is analytic throughout the real plane. Then 
a necessary condition for this function to take a minimum 
is d^f - dyf - 0. It is easy to check that dyf is a sim- 
ple second-degree polynomial in //, with zeros ji^i - Q and 

^^.2^iJ\+yl-^7l+^xl^)|[iy]-yl+xl^)+iyl-y\)vl The former 
solution clearly cannot correspond to stationary points of / - 
in particular zeros of/, i.e., singular points ofh - that lie in the 
interior of ^S (as anticipated, a zero of/ is also a minimum and 
thus one of its stationary points). On the other hand, by plug- 
ging n^2 into <9^/ we find df,f\ft=^^^ ^'^i7^i-^l+xlj,)iyl-yl+xl^). 



which vanishes for either y^-jl+x^^ -^oi y\-y\+x\-^ - 0. We 
have already shown (see above) that in neither of these two 
cases h can admit minima in the interior of S. 



B. Minima on the boundary of S 

The findings of the previous subsection show that we can 
restrict the search for the minimum of h to the boundary of S. 
The possible values of h on the square edges corresponding to 
ji-Q, ji-l, v-Q and v = 1 are, respectively, given by 



K=o^y\+yl+\y\-yl\^^y\, 



(72) 



h^=i ^yl+x\^+y\+{j\-y\)v+\yl-yl+xl^+{y\-yl)v\ , (73) 
hy=Q^yi+y\-iy\-yl-x\^)^i+\yl-y\{l-^i)+(xl^-yi)^i\ilA) 
hy=l^yl+y\-iyl-yi-x\^)^i^y]-yl{l-^i)+{x\^-yl)^i\.{15) 



From Eq. ( [72] i it trivially follows that the minimum of h on 
edge // = is given by min /z^=o - ^y^. In the next three dedi- 
cated paragraphs, we minimize h on edges fi-l and v = 0, 1, 
respectively. 



1. Edge yu = 1 



This is the set of points (/i=l,0<v<l)on which function 



h is given by Eq. (73 1. Let h+ (hJ) be the expression taken by 



h when the absolute value in Eq. (73 1 is positive (negative). 
These are easily calculated as 

hAv)^2[yl+x\, + (y\-yl)v\, hjv)^2yl. (76) 

Importantly, note that h+ always grows with v while /!_ is flat. 
The argument of the absolute value [cf. Eq. ( [73) l] increases 
with V (since y^ > yf) and vanishes for v = vq = ~iy\~y\ + 
x^3)/(7j -72). Hence, it is negative (non-negative) for v < vo 
(v > Vo). As a consequence, h-h^ (h = h+) for v < vo (v > vq). 
Now, the minimum of h on this edge depends on the sign of 
Vo, which depends in turn on the sign of yj-yj-i-x^j. Indeed, 
if 72 - 73 -t- x^ J < then vo > and thus min /i;,= i = min /;_ = 27^ 
(recall that h+ grows with v). If, instead, y2~yl+^A3 - ^ 
then Vo < and h = h+ for < v < 1, namely throughout the 
edge. The minimum is thus taken at v = and reads min/i^=i = 
min/!+(v = 0) = 2(72-i-x^3). To summarize. 



if yl-yl+xl3<0 min/i^=i=2y^. 



(77) 



if 72-r3+-^l3^0 min /i^=i =2(7^-^x^3), (78) 



2. Edge V = 



This is the set of points (0</i< 1, v = 0), where /1 is given by 



Eq. (74 1. Similarly to the previous paragraph, we first search 
for the zero of the absolute value, which is easily found as 
ji- jiQ- iy]-yf)/(y^-yl+x^^). Its location on the real axis 
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fulfills 



ifrf 



-r^ 



■^A3 



>0 



ifr?-r3+-^A3<0' 



ii7l-jl+xlj>0, 0<jUo<l, 



if 72 
;U0<0 , 



73+^A3<0' 



A*!) > 1 . 



(79) 
(80) 



yields h(fj. < fio) = /;_ and h(fi > /iq) = h+. Now, two cases 
can occur. If \yj,\ > \yi\, then h- grows with yu and therefore 
min /Zv=o - h-(fJ. = 0) = 2-/^. If |y3| < lyj |, instead, /z_ decreases 
with yU. Then the minimum of h depends on whether or not 
fio^ 1, which depends in turn on the sign of y^-y^+x^^ ac- 



cording to Eq. (79 1. If 72~ yi'^^A^ - then /yq < 1 and h is 



minimized for// =/i() (recall that h+ always grows). This yields 



which we will use in our analysis. At variance with the previ- mmh,.=o = /i-(A'o) = 2 [yf(y^ + x^3)-r^r^]/(7f-r^+x^3). On 

the other hand, j\-y\-^ ^/^-t, < implies jUq > 1. Hence, 
h = /i_ throughout the interval < // < 1 and, necessarily. 



ous paragraph, now the absolute value in Eq. (74 1 grows (de 



creases) with/i ioiyf-y^+x'^^>Q (yf-y^+xj^^KQ). Eq. (73 i 
straightforwardly gives 



min hy^o = h_{ji-\)- ly^. 



Kifl) = 2{yl + xl,^i\ /2» = 2 [y] + {yl-y])^j\ , (81) 

where h±- are defined in full analogy with the previous para- 
graph. Note that, while h+ always grows with /u, h^ is an 
increasing (decreasing) function of /( for |y3| > \y\ \ (|y3| < lyi |). 
Let us analyze the possible situations. Based on the above, 
if -/j - y^ -I- jc^3 > then jUq > and, moreover, the absolute 
value is negative (non-negative) for fi < i-io (fi > juo)- This 



< 0. In this situation. 



We are left with the case J^-yl+x^^ 
fi[) < [cf. Eq. ( [80| ] and the absolute value is non-negative 
(negative) for ji < jio {p.> /zq), which gives h = h^ through- 
out this edge. Now, the analysis is simpler since, evidently, 
only the case \yi\ < \yi\ is possible. Thus /;_ can only increase 



[recall Eq. (|8 Ik] and min hy^o -h^iii-Q)- 2y^. 



To summarize, on the edge v = 



ifr?-r3+43<o 



ify\-yl+x\^>0 



mmhy^Q-2y^ , 
' if 1731 > In I minhy^o^2y^i , 

if|73l<l7il min/iv=o = 0(72-73+-^L) 2 



7 1 "73 +^A3 



+ 0[-(72'-7^^L)l27 



(82) 
(83) 



3. Edge V = 1 



This is the set of points (0</i<l,v=l), where /i is given by 



Eq. ( |75) . Similarly to the previous paragraph, we first search 
ero of the ab 

-(7?-7^)/(7^ 



for the zero of the absolute value, which is easily found as 

jj-fj^o = -(7i -72)/(72~73+-^A3)- ^'■^ location on the real axis 
fulfills 



if 72-73+^3^0 



//o<0. 



(84) 



if yl-yWA,<ol''/nt'h^l 0<^f 1 (85) 

n /3 A3 ^jf^2_^2+^2^<Q ^^^>^ 

Based on Eq. ( [75| , the expressions taken by h on this edge 
when the absolute value is positive and negative are, respec- 
tively 



7^ 



yl+x^2 <0 then 



thus min/iv=i - h+{Q) = 2y'^. If, instead, 
h-h+ (h-h-) for fi<fi() (/OyUo). Moreover, note that in this 
case one has y^ > y^, which entails [cf. Eq. ( 86 1] that both /;_ 
and h+ grow with //. Hence, the minimum is again given by 
min hy^i =h+{Q) = 2y^, which completes our proof. 



C. Global minimum 

To give the general expression for the minimum of h it is 
convenient to refer to the minimization study on the edge v = 0. 
Recall that the minimum of h on the edges fi-Q and v = 1 is 



unconditionally given by 2yj. When 77-73 +-^^3 < 0, based 
on Eqs. ( 77 1 and ( 82 1 the minimum reads min^ h = 2y^ (note 
indeed that this case necessarily entails y2-y\ + x\^ < and 



K{^i) = 2iy\+xl^^i), /2» = 2 k' + (73-72)j"l ■ (86) rf <73)- If instead y\-y'^+x\^>0, both signs of yf-y^ 



-,+x 



A3 



Hence, h+ always grows with ^ while h^ is an increasing (de- 
creasing) function of // for 173! > |y2l (I73I < l72l)- We show next 
that, the minimum of h on this edge is always given by 2y\. 



as well as 1731-1711 are possible. Hence, if I73I > I71I upon 
analysis of Eqs. ( [77] l, ( [78| and the first case in Eq. ([83| ) we 
end up with ming h - min{27p 2{y^+x\^)). For T^-J^+^^t, < 



0, this gives min^ h 



Indeed, if 72^73 +-*^a3 - ^ (implying y\- 
fio < and h = h+(fj.) throughout the edge. 



73-i-x^3 >0) then 
The minimum is 



27^ 

St 

2[7f(7^43)-7^73']/(7^ 



For y\-y\ + x^ 



A3 



> 0, the global 



minimum is the lowest number among 2y^^, 2{y2+x^^) and 



yl+x 



l,y 



Hence, to summarize. 
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'■if|73l>lril mmh^&{yl-yl+xlj)2mm{ylyl+xl^}+&[-{yl-yl+xl^)]2yl. 



ifyi-yi+xi^>0 



iflyaKlnl min /z ^ Q {yj -y^xi3) 2 min (y^ yl+ x^^ '^^'^T^'^ }+0 [- {yl -yj +xl,)] lyj. 



Eq. ( [87|) can be further simplified. Indeed, on the second line 

(case^-y^ + ^c^j > and lyjl > |yi|) for yl-yj + xl^ > 
we have y^+x^^ >y\>y\ and therefore the minimum is 2y\ 

On the other hand, on the third 



713 
.2 
-A3 



^r^ 



regardless of yj + . 
line (case y\-y\ + - 

and using y\ > y\ it is straightforward to prove that the ratio- 
nal function can never exceed both y\ and y\ + y?}^^- In light 



>0 and lysl < |yi|) for y2-r3+-^A3 



of these considerations and upon comparison of Eq. ( 87 1 with 



Eqs. (82 1 and (83 i, we conclude that the global minimum of 
h is achieved on the edge v = 0. Therefore, using Eq. ( 32 1 



the TDD of an arbitrary two-qubit X state is finally obtained 



as in Eq. ( 65 1. Remarkably, in each case that can occur (de- 
pending on the parameters defining the state) 'D^^\pab) takes 
a relatively compact expression. 

As already anticipated, for Bell-diagonal states (see Section 



nil, Eq. (65 I yields the result of Section 111 A as shown in 
detail in Appendix B. 

Another interesting special case occurs when in Eq. ([63]l 



either p32 = or p4i = 0. Then, due to Eq. (64i, |yi| = \y2\ 
-^^ — — — ^-^ ^ ^. Hence, such a case always entails 



and 



D^^^(Pab) - ItiI/2, namely half of the absolute value of the 
non-zero off-diagonal entry. 



VI. CONCLUSIONS 

In this paper, we have addressed the issue of the com- 
putability of TDD, one of the most reliable and advantageous 
QCs indicators. By introducing a new method for tackling 
and simplifying the minimization required for its calculation 
in the two-qubit case, we have demonstrated that this can be 
reduced to the search for the minimum of an explicit two- 
variable function. Then, we have shown that this can be an- 
alytically found in a closed form for some relevant classes 
of states, which encompass arbitrary quantum-classical and X 
states. The latter includes as a special subset the Bell diago- 
nal states, which were the only states for which an analytical 
expression of TDD had been worked out prior to our work. 

Due to the importance of quantum-classical and X states, 
along with the typical hindrances to the calculation of their 
QCs through bona fide measures, our work provides a signifi- 
cant contribution to the study of QCs quantifiers, by combin- 
ing the desirable mathematical properties of TDD with its ex- 
plicit computation for these classes of density matrices. Fur- 
thermore, we expect that the framework developed in this pa- 



(87) 



r 



per may be further exploited in future investigations to enlarge 
the class of quantum states that admit an analytical expression 
for TDD. 
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Appendix A: Derivation of Eq. J48| 

We recall that f±i(p, a) -y cos (f> + xa cos,((f> - a). This is a 
linear combination of cos (p and sin (p, which can be arranged 
in terms of a single cosine as A-t(cos (p cos 6± + sin (p sin 6±) = 
A-t cos{<p-6±). Using xai -xa cos a [see Eq. (44 1], the factor 

is easily found as A± - Jiy+XA\)^+x^2^ while 



sin^-t 



cos6± 



±xa sm a 



±Xa2 



^{y + xax)^ + x\2 ^(y + xai)^ + 



(Al) 



2 _J_ y2 

^A2 



y + Xa cos a 



y±XAi 



(A2) 



^^-^L 



^(y + xax)^ + x\2 ^{y±XA\9 + 
Hence, (5±=arctan[xA2/(.JAi±T)]- Therefore 

f±{^, a) = y (y + XA\9- + x\2 cos(cf) - (5±). (A3) 



Replacing Eq. ( |A3| l into Eq. (47 1 of the main text yields 
\\M(e)\\i^2j] yl[(y+ijXAi)^+xl2] [l-cosHcf>-5,)] 

^2Y^yj(y+i]XAi)^+xl2 \sm(4)-6,;)\ 

-2y ■J(y+riXAi)^+x^2 I sin cos (5;, -cos sin 5,, 



i=± 



Eliminating now sin 6-t and cos 6-t through Eqs. (All and ( A2 1 



we end up with Eq. (48 i of the main text. 
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Appendix B: Eq. J65) for Bell diagonal states 



172! and the TDD reduces to 



Bell diagonal states are defined as a mixture of the four 
Bell states. This immediately yields that they fulfill xa-xb- 
0, that is, the reduced density matrix describing the state of 
either party is maximally mixed. Therefore, the corresponding 
density matrix can be expanded as a linear combination of 
lyi®!^ and {a-Ak®crBk]- As each of these four operators has 



an X-form matrix representation [cf. Eq. ( 63 1] Bell-diagonal 
states are X states. Hence, in Eq. (65 1 ■yj-YY^x^^^yf-y^ for 
i-\,2. In this case, the square root m Eq. ([65]l coincides with 



Ir3l>lril ^ £'*^Vab) = 



In 



1 



IraKlnl ^ 0'"\pab)^ :^raiix{\y2\,\yi\ 



(Bl) 



It is immediate to check that the above is equivalent to 
state that D^'^Xpab) is half of the intermediate value among 
{ItiL IT2L ItsI), which fully agrees with Refs. llT2l[T3l and the 
findings of Sec. [Ill A| 
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